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Abstract 

We consider discontinuous dynamics of integrate-and-fire models, which consist of pulse-coupled biological 
oscillators. A thoroughly constructed map is in the basis of the analysis. Synchronization of non-identical 
oscillators is investigated. Significant advances for the solution of second Peskin's conjecture have been 
made. Examples with numerical simulations are given to validate the theoretical results. Perspectives are 
discussed. 
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1. Introduction 

The collective behavior of biological and chemical oscillators is a fascinating topic that has attracted a 
lot of attention in the last 50 years |lj]-[28j]. An exceptional place in the analysis belongs to synchronization, 
which in its general sense is understood as phase locking, frequency locking, and synchrony itself, that is 
motion in unison 0, Q,@"[23- The integrate-and-fire model of the cardiac pacemaker 34 1 was developed 



by C. Peskin 35] to a population of identical pulse-coupled oscillators. It was conjectured that the model 
self-synchronizes such that: 

(CI) For arbitrary initial conditions, the system approaches a state in which all the oscillators are firing 
synchronously. 

(C2) This remains true even when the oscillators are not quite identical. 

The conjecture (CI) is solved in 35] for a system with two oscillators, and in [2l[ for the generalized model 
of two and more oscillators. The last paper gave start to an intensive and productive investigation of the 
problem and its applications 0,1, El, [36|-[4i[. As far as we know the conjecture (C2) remains unsolved. 



Even a developed non-identity concept has not been found in the literature. 

In the present paper we generalize the model, and propose a version of non-identity. The model is 
considered such that perturbations save the synchronization. These oscillators are not only pulse-coupled, 
but connected during the time between moments of firing. That is, the modeling differential equations are 
not separated. One can see that this approach may provide more biological sense to this theory. The paper 
consists of main results, simulations and discussion of possible extensions. We believe that the results and 
proposals of this paper reveal new perspectives on the study of integrate-and-fire models of oscillators. The 
main role in our analysis is played by a specially defined map. It is not a Poincare map, since it transforms 
the coordinate of one oscillator to that of another, and the two interchange roles in the course of the mapping. 
If there are more than two oscillators, they are used in pairs to shape the map with the interference of other 
oscillators acting as perturbation. 
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The main object of our investigation is an integrate-and-fire model, which consists of n non-identical 
pulse-coupled oscillators, Xi,i — 1,2, ... ,n. If the system does not fire the oscillators satisfy the following 
equations 

x[ = f( Xl ) + Mx)- (1-1) 

The domain consists of all points x — [x\, X2, ... , x n ) such that < xi < 1 + Ci( x ) f° r au * = 1,2, ... ,n. 
When the oscillator Xj increases from zero, and meets the surface such that Xj(t) = 1 + Q(x(t)), then it 
fires, Xj(t+) = 0. This firing changes the values of all oscillators with i ^ j, 

Xl[t+> - { Xi(t) + e + ei, otherwise. [ > 

Thus, it is assumed that if Xi(i) > l + Q(x) — e — €i, then the oscillator fires, too. It is assumed also that there 
exist positive constants fii and ^ such that < fa and |Ci(^)| < £ij f° r au x an d I n what follows, we 

call the real numbers e, fa,£i,£i, parameters, assuming the first one is positive. Moreover, constants £j, e,, fa, 
will be called parameters of perturbation. If all of them are zeros, then the model of identical oscillators is 
obtained. We assume that e + e-i > 0. That is, an exhibitory model is under discussion. The function / is 
positive valued and lipschitzian. Moreover, assume that Q are continuous and fa are locally lipschitzian for 
all i. 

The coupling in the model is all-to-all such that each firing elicits jumps in all non-firing oscillators. 
If several oscillators fire simultaneously, then other oscillators react as if just one oscillator fires. In other 
words, any firing acts only as a signal which abruptly provokes a change of state. The intensity of the signal 
is not important, and pulse strengths are not additive. A system of oscillators is synchronized if all of them 
fire in unison. 

In the present analysis we address synchronization as well as the existence of periodic solutions. Results 
that concern continuous and delayed couplings are considered in [42, 43]. 

We believe that the approach proposed in this paper will be useful for the investigation of a wide range 
of problems, focusing not only on synchrony and pulse-couplings, but also phase locking, frequency locking 
of systems, families of oscillators with continuous couplings. The method can be used to analyze inhibitory 
models as well as to evaluate the effects of coupling time deviations. Moreover, the model is suitable for the 
investigation of the existence of quasi-periodic and almost periodic motions. 

2. The prototype map and two identical oscillators 

In this section we shall define the map, which is the basic instrument of our investigation. It is constructed 
for a model more general, than is needed for this paper, to be the basis for future investigations. 

Let us consider two identical oscillators, x\{t),X2{t),t > 0, which satisfy the following differential equa- 
tions 

A = f{xi), (2.3) 

where < < l,i = 1,2. When the oscillator Xj fires at the moment t such that Xj(t) — l,Xj(t+) = 0, 
then the value of another oscillator with i ^ j, changes so that 

Xl = { x\{t) +e, ''otherwise. 
Denote by u(t, 0, uq), the solution of the equation 

u = f(u), (2.5) 

such that u(0,0,i4o) = Uo- Assume that the solution exists, is unique and continuable to the threshold for 
all uo- Consider the solution u(t) — u(t,0,v + e) of (|2.5I) . Denote by s(v) the moment when u(s) — 1, and 
define the function L(v) — u(s, 0, 0) on (0, 1 — e). 

The following conditions will be needed throughout the paper: 
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(Al) L(v) is a strictly decreasing continuous function; 
(A2) t] = lim^ 0+ L[v) > 1 - e; 
(A3) linw_ e L(v) = 0. 

Conditions (Al), (A3) are valid, if, for example, / is a positive and lipschitzian function. Another case will 
be considered in Example 12.31 It is obvious that there exists a unique fixed point, v*, L(v*) = v* . 
Now, define a map L : [0, 1] [0, 1], such that 



This newly defined function is continuous on [0, 1]. The sketch of its graph is shown in Figure [TJ To make 



Figure 1: The graph of function w = L(v), in red, and the period— 2 orbit in blue. The points ao = 0,ai = 1 — e,afc +1 = 
L~ 1 (a k ), k = 1,2, 3. (Color online) 

the following discussion constructive consider the sequence of maps L k (v),k — 1,2,..., where L k (v) = 
L(L k ~ 1 (v)) if k > 2. Their graphs with k — 1,2,3 are shown in Figure^] Denote ao = 0, ai = 1 — e, a2 = 

L _1 (l — e),a 3 = (L 2 ) _1 (l — e), The sequence can be obtained also through iterations a = 0,ai = 

1 — e, afc+i = L^ 1 (ak), k = 1,2,..., which are seen in Figure[TJ It is clear that the sequences an and ci2i+i 
are monotonic, increasing and decreasing respectively. One can verify existence of a fixed point v** < v* 
of the map L 2 (v) such that v = L(v**) > v*, and there are no fixed points of L 2 in (0,v**). Moreover, 
Q-2i ~^ v ** an d a2i+i - 1 v as i — > oo. In the case i>** = v* , there is no non-trivial period— 2 points of L. 

Let us show how iterations of L can be useful for the investigation of synchronization. Consider a motion 
(xi(t),X2(t)) and a firing moment t$ > such that x\(to) = l,xi(to+) = 0,^2(^0+) = v,v G [0, 1]. 

Lemma 2.1. Motion (x\(t),X2(t)) synchronizes if and only if there exists a number k such that 1 — e < 
L k (v) < 1. 

Proof. Let us consider only necessity, since sufficiency is obvious. We shall consider the following two cases: 
(a) < v < 1 - e; (f3) 1 - e < v < 1. 

(a). It is clear that the couple does not synchronize at the moment t = to. While it is not in synchrony, 
there exists a sequence to < t\ < . . . , such that x\ fires at moments ti with even i and X2 fires at i$ with 




L(v), if we (0,1-e), 

7/, if v = 0, 

0, if [l-e,l]. 



(2.6) 




3 



Figure 2: The graphs of L, L 2 and L 3 in red, blue and green respectively. (Color online) 



odd indexes. Set Vi = xi(U) if i is odd, and Vi = £C2 (ii) if it is even. Use the definition of L and identity 
of oscillators to obtain that Vi + \ — L(vi),i > 1. This demonstrates that map L evaluates alternatively the 
sequence of values X\ and xi at firing moments. 

The pair synchronizes eventually if and only if there exists k > 1 such that x\(t) ^ Xi(f), if t < tk, and 
xi(t) = X2(t), for t > tk- Both oscillators have to fire at tk- That is, 1 — e < Vk < 1. 

(j3). Consider 1 — e < v < 1, as the case v = 1 is primitive. We have that to is a common firing moment 
of both x\ and x%, and it is the synchronization moment. Moreover, 1 — e < L 2 (x2{t\)) — r\ < 1. The lemma 
is proved. □ 

Thus, it is confirmed that the analysis of synchronization is fully consistent with the dynamics of the 
introduced map L(v) on [0, 1], and, therefore, L can be used as a valuable tool in further investigation of 
the topic. 

Let us consider the rate of synchronization. We solve the problem by indicating initial points which 
synchronize after precisely k, k > 0, iterations of the map. Denote by Sk the region in [0, 1], where points v are 
synchronized after k iterations of map L. One can see that So = [1 — e, 1], S\ = [ao, 02], and Sk — (dfc-i, £ifc+i], 
if k > 3, is an odd positive integer, and Sk = [flfe+i, dk— l); if Jb > 2, is an even positive integer. 

From the discussion made above it follows that the closer v is to v** from the left or to v from the right, 
the later is the moment of synchronization. 

Denote by T the natural period of oscillators, that is, the period, when there are no couplings, and by 
T the time needed for solution u(t, 0, v*) of (I2.5[) to achieve threshold. Since each oscillator necessarily fires 
within any interval of length T and the distance between two firing moments of an oscillator are not less 
than T, on the basis of the above discussion, the following assertion is valid. 

Theorem 2.1. Assume that conditions (Al) — (A3) are valid, and to > is a firing moment such that 
xi{to) = l,Xi(to+) = 0. If X2(to+) G S m ,m is a natural number, then the couple x\,X2 synchronizes within 
the time interval [to + ', to + mT\. 

One can easily see that whenever condition (A2) is not valid, the system does not synchronize. 

Example 2.1. Consider the integrate- and- fire model of two identical oscillators Xx,X2, with the differential 
equations 
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(2.7) 



where i = 1, 2, and c is a positive constant. It is known \9^] that the canonical type I phase model Jj/ can be 
reduced by a transformation to the form 

u' = u 2 +c. (2.8) 

This time we investigate the model with the pulse- coupling. 

Since the two equations are identical, we shall consider a solution u(t) of equation \2. 8\) to construct map 

L. We have that u(t, 0, v + e) = -y/ctan(ci + arctan( !i iS)) and 

v c 

V ~\~ 6 

V^tanfcs + arctanf — ^)) = 1- (2-9) 

Next, u(s, 0,0) = •y / ctan(cs), and by applying 12. 9\) we find that 

1 — v — e 



L(w) = c- 



if v E (0, 1 — e), and the fixed point is equal to v* — y/(c + e/2) 2 + c(l — e) — (c + e/2). Evaluate 

1 - e 



L(0) 



c + e 



<o see that L(0) < 1 — e, and condition (A2) is not valid. Moreover, one can verify that L'(v) < 0. 
Thus, we obtain that the couple does not synchronize, and our simulations confirm this. 

Example 2.2. Consider the following integrate- and- fire model of two identical oscillators, X\^x%, with the 
differential equations 

x'^S-^Xi, (2.10) 

where i = 1,2, positive constants 5,7 satisfy k = — > 1. One can find that u(t,0,v + e) = (v + e)c ft + 
k(1 — e~ 7 *) and u(s, 0, 0) = k(1 — e _7S ). T/ie Zasi <wo expressions imply that 

L(,) = J^±4, (2.11) 
w K - (w + e) v ' 



ifO < v < 1 - e. 

There is a unique fixed point of L and L 2 , and it is equal to 



„* = («--)- A /«2_, c+ _ (2.12) 



Finally, L(0) = k-^— f > 1 — e. T/iai is, a/Z conditions of the last theorem are valid, and the assertion in \3, 
is proved. 

Remark 2.1. Map L is similar to that in fSa] . but the argument here is a coordinate before a jump, while 
in the paper the argument is a coordinate after a jump. This difference is not a critical one. The most 
important point is that C. Peskin uses it only as an auxiliary device to build the Poincare map. We use L 
itself as the main map, with a newly defined continuous extension, which simplifies the discussion in this 
section and throughout our investigation. We revisit the problem of two identical oscillators, since L is the 
prototype map in our analysis. In addition to the main synchronization result, regions with equal time of 
synchronization are indicated, and the value of the fixed point, v* , is evaluated. 
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Example 2.3. Consider the following integrate- and- fire system of pulse- coupled oscillators, Xi,x 2l such that 

x 'i = f(xi), 

where 

( 4 -3s, if < s < 1/3; 

f(s) = \ 3, if 1/3 < a < 2/3; 

[ 4- 3(s- 2/3), if 2/3 < s < 1. 

We have found that map L for this system exists and is equal to 

r 2§±|f, if o<^<i/3- e; 

L(v) = < 1-v-e, if 1/3- e < t> < 2/3- e; 
I ffcf, if 2/3- e <«<l. 

One can check that conditions (Al) — (A3) are fulfilled for this map. Moreover, the fixed points are equal 
to v* = (1 — e)/2,v** — 1/3 and w = 2/3 — e. Finally, all the motions, which start outside of the periodic 
trajectory synchronize eventually, and all of them are periodic inside the trajectory. 

The last example shows that the assumptions for map L, including the existence of the non-trivial period— 2 
motion, can be realized for even the differential equations with discontinuous right-hand side. Moreover, 
in future investigations one can consider isolated periodic solutions, stable or unstable. The theoretical 
consequences of this research are clear, if one uses the mappings theory, but construction of examples 
requires additional time. 

Example 2.4. Consider the model of two integrate- and- fire identical oscillators, x\,x 2 , which are pulse- 
coupled and 

x[ = S - 7x1 + (3x 2 , 

x' 2 = S -jx 2 + /3 Xl , (2.16) 

where constants S, 7 and f3 are positive numbers. One can easily see that the system is the extended Peskin's 
model in Example \2.2l The terms with coefficient j3 are newly introduced in the system. They reflect the 
permanent influence of the partners during the process. Eigenvalues associated to \2.1b}) are \\ = —7 + j3 
and \ 2 = —7 — f3. We suppose that f3 is small so that both eigenvalues are negative. Moreover, it is assumed 
K = 5/7 > 1. Then, K\ = —S/\\ > 1 if (3 is sufficiently small. The solution of system \2. 1 6}) with value 
(0,«+e) att = 0, is equal to ui(t) = 5 [e Al '-e A2t ](v+e)-Ki(e Alt -l), u 2 (t) = ±[e Alt +e A2t ](v+e)-Ki(e Alt -l). 
By using these expressions obtain the equation i[e AlS + e X2S ](v + e) + Ki(l — e Als ) = 1, and construct 
L(v) = 1 — (v + e)e A2S . Map L is too complex to analyze for properties (Al) — (A3). That is why we will 
compare this model with the couple in Examvle \2.2\ The last two equations imply L(Q) = K^rf > 1 — e, if 
(0 = and v = 0. That is, if f3 is sufficiently small, then condition (A2) is valid. We have found, also, by 
direct evaluations that the derivative L'(v) is negative if S and [3 are sufficiently large and small respectively. 
That is, condition (Al) is fulfilled. It is easy to verify that condition (A3) is also correct. 

Now, using the continuity theorem in parameters fJ2] . one can find that map L may admit a period— 2 
point only if the orbit is as close to the fixed point v* of Example \2.2\ as (3 is small. Consequently, the 
measure of the set of points, which can not be synchronized diminishes as (3 — > 0. This result is a new one. 
In previous papers the differential equations were separated. 

3. The multidimensional system of non-identical oscillators. 

Consider the model of n non-identical oscillators given by relations (|1.1[) and (|1.2[) . The domain of this 
model consists of points x — (x\,x 2 , ■ ■ ■ , x n ) such that < Xi < 1 + Q(x) for alH = 1, 2, . . . , n. 
Fix two of the considered oscillators, let say, x\ and x r . 
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(2.13) 
(2.14) 

(2.15) 



Lemma 3.1. Assume that condition (A2) is valid, and to > is a firing moment such that xi(to) = 
1 + Q{x(to)), Xi(to+) — 0. If parameters are sufficiently close to zero, and absolute values of parameters of 
perturbation sufficiently small with respect to e, then the couple xi,x r synchronizes within the time interval 
\t$,tQ+T\ ifx r (to+) [ao,di) and within the time interval [fo + " 1 ^" 1 T, to + (m+l)T], ifx r (to+) G S m ,m > 1. 

Proof. Denote by x(t) — [x\(t) , x 2 {t) , . ■ ■ ,x n (t)), the motion of the oscillator. If 1 + ( r {x(to)) — e — 
e r < x r (to) < I + Cr(x(to)), then these two oscillators fire simultaneously, and we only need to prove the 
persistence of synchrony which will be done later. So, fix another oscillator x r (t) such that < x r (to) < 

1 + ( r (x(t )) -£-S r . 

While the pair does not synchronize, there exists a sequence of moments < to < ti < • • ■ > such that 
oscillator xi fires at ti with even i, and x r fires at ti with odd i. For the sake of brevity let Ui = xi(ti),i = 
2j + l,Ui = x r (ti), i = 2j, j > 0. In what follows we shall evaluate the difference itj+i — L(iii). 

Let us fix an even i and Ui = x r (ti). If the parameters are sufficiently small, then there are k < n ~ 2 
distinct firing moments of the motion x(t) on the interval (ti, U + \). Denote by U < 0\ < 02 < ■ ■ ■ < 0^ < U + \, 
the moments of firing, when at least one of the coordinates of x(f) fires, and v(t,to,Vo) the solution of the 
equation (jl.ll) with v(to,to,vo) =Vq. We have that 

x r (9i) = x r (ti) + e + J f(x r (s))ds + 



(j) r (x(s))ds, (3.17) 



where x(t) = v(t,U, x(U+)), 



x r (9 2 ) = x r (6i) + e + / f(x r (s))ds + 

<f> r (x(s))ds, (3.18) 



where x(t) — v(t,9i,x(9 1 +j), 



rU+i 

x r (U +1 )=x r (6 k ) + €+ / f(x r (s))ds + 

/ <t> r (x(s))ds, (3.19) 

where x(t) = v(t,6\,x(6k+)). 
The moment ti+i satisfies 

1 + Cr(ar(ti+i)) -£-£,-< ar r (*<+i) < 1 + Cr(*(*i+i)). (3.20) 
Similarly to the expressions for x r one can obtain 
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xi{6x)= j f(x l (s))ds+ [ 1 </n{x(s))d8, 
Jti Jti 



x l (9 2 )=xi(9 1 ) + e+ I f(x l {s))ds + 
4>l(x(s))ds, 



fU+i 

xi{t i+ i) = xi(9 k ) + e + I f{xi(s))ds + 



<t>i(x(a))ds. (3.21) 

Formulas (|3.17l) to (|3.21[) define Mj+i = xi(ti + i). Similarly one can evaluate the number for odd i. 
Let us now find the value of L(ui). With this aim, evaluate 



<P(i i+ i) = x r {U) + e + / f(4>(s))ds, (3.22) 

where £j + i satisfies = 1, and 

m+i) = I f(ip(s))ds, (3.23) 



to find that L(iii) = "0(^+1 )• Next, we will show that the difference u 2 +i — L(ui) is small if the parameters 
are small. 

First, one can find 

4>(t)-x r (t) = 

[f(4>(s))- f(x r (s))]ds- f Mx(a))da, (3.24) 



for t G [U,0i]. 

Then, by applying the Gronwall-Bellman Lemma one can easily see 

10(0!) - av(ffi)| < - tije^ 1 "* 4 ), (3.25) 
where I is the Lipschitz constant of /. Next, we have 

\<f>{9 2 ) - x r {9 2 )\ < 

\Hr{Bx - U)e e ^-^ + fi r (9 2 - Ox) + ( ]S ^\ (3.26) 

if te [9 1 ,e 2 \. 

Without loss of generality, assume that t i+1 > ti + \. Proceeding the evaluations made above, we can 
obtain |1 - i r (t i+1 )| = \4>{t i+1 ) - x r (i i+ i)\ = $i(e, /x r ), where 

fc-i 



$i(e,iM.) = M(#i - U)S 1 ^-^ + ^(^+i - 9j)e e ^- 



3=1 

(t i+ i - fffc)^^ 1 - ^] + e ^e^+ 1 -^). 

3=1 

There are positive numbers fi and M, which satisfy /i < f(s) < M, if < s < 1+max^ £j. One can request the 
following inequality maxj = i i ... i „ fii < fi. We have that \x r (ti + i) — x r (U + i)\ < |1 — x r (tj+i)| + |1 — x r (ii+x)\ < 
$i(e,/i r ) + £ r . Consequently, 

$i(e, u r ) + £ r . ^ . 



By applying (|3.2ip and (|3.23[) . making similar evaluations for (|3 . 25[) and (|3.26p one can find \ip(ti+i) — 
xi(ii+i)\ < $i(e,w)- 

Then, we have that \u i+1 - = \ip{t i+ i) - xi{t l+1 )\ < \ip(t i+1 ) - xi{t i+1 )\ + \xi(t i+ x) - Xi(t i+1 )\, 

and, consequently, 

- L( Ui )\ < $(e, fi r , w , &.), (3.27) 

where $ = $1 + + [i r ). It is obvious that $ tends to zero as the parameters do. This convergence 

is uniform with respect to uq. We can also vary the number of points 9i and their location in the intervals 
(tj,tj+i) between and n— 1. The convergence also is indifferent with respect to these variations. 
Consider the sequence of inequalities 

|«i - £>o)| < \ui - + \L{ui- X ) - L(L i - 1 (u Q ))\,i= 1,2, ... . 

Then recurrently, by applying continuity of L, (|3.27[) and L rn (uo) G [1 — e, 1], conclude that either 1 + £ r — 
e — e r < u m < 1 + £ r or 1 + — e — e; < u m +i < 1 + 6: if the parameters are sufficiently close to zero, 
and absolute values of the parameters of perturbation are sufficiently small with respect to e. Both of these 
inequalities bring the pair to synchronization. 

Since each of the iterations of map L happens within an interval with length not more than T, we obtain 
that couple xi, x r synchronizes no later than t = to + (m + 1)T. Similarly, the couple synchronizes not earlier 
than t = t + ^f. 

If two oscillators x\ and x r are non-identical and fire simultaneously at a moment t = 9, how will they 
retain the state of firing in unison, despite being different? To find the required conditions, let us denote 
by t,t > 6 a, moment when one of them, let's say x r , fires. We have that xi(6+) = x r (8+) — 0. Then 
Xi{t) = x r (t),9 < t < t. It is clear that to satisfy xi(t+) = x r {r+) = 0, we need 1 + C,i{x{t)) — e — e; < Xi(t). 
By applying formula (|3.20[) again, this time with ti = 9, ti+\ — r, one can easily obtain that the inequality is 
correct if parameters are close to zero, and absolute values of the parameters of perturbation are sufficiently 
small with respect to e. Thus, one can conclude that if a couple of oscillators is synchronized at some moment 
of time then it persistently continues to fire in unison. The lemma is proved. □ 

Remark 3.1. The last lemma not only plays an auxiliary role for next main theorem, but can also be 
considered a synchronization result for the model of two non-identical oscillators. 

Let us extend the result of the lemma for the whole ensemble. 

Theorem 3.1. Assume that condition (A2) is valid, and to > is a firing moment such that Xj(to) = 
1 + Q(x(to)), Xj(to+) = 0. If the parameters are sufficiently close to zero, and absolute values of parameters 
of perturbation are sufficiently small with respect to e, then the motion x(t) of the system synchronizes within 
the time interval [to,to+T], if Xi(to+) $ [ao,ai),i ^ j, and within the time interval [to + max '^ 2 J kl — -T,to + 
(max^j ki + 1)T], if there exist a; s (£o+) € [do, o-i) for some s ^ j and Xi(to+) € ,i ^ j. 

Proof. Consider the non-trivial case. Applying the last lemma we can see that each pair (xj,Xi),i ^ j, 
synchronizes within [to + max '?^ fc » — -T,to + (max^j ki + 1)T"]. The theorem is proved. □ 
Now, replace coupling (|1.2[) by 

*(*+) = { 0. ^(*)+* + ^l + 6(»). (3.28) 
v ' [ Xi{t)+e + ei, otherwise, v ' 

where e, e + ej > 0, is a new parameter, independent of e. Consider system (jl.ip with p.28p . One can find 
that the following assertion is valid. 

Theorem 3.2. Assume that condition (A2) is valid, and to > is a firing moment such that Xj(to) = 
1 + Q(x(to)), Xj(to+) = 0. If parameters e,^i,^,ej, are sufficiently close to zero, then the motion x(t) of 
the system synchronizes within the time interval [to, to + T], if Xi(ta+) [aa,ai),i ^= j, and within the 
time interval [to + max '^ fc ' — -T, to + (max^j ki + 1)T], if there exist x s (to+) G [ao, ai) for some s ^ j and 
Xi{t +) G Sk it i ^ j. 
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We can see that (|3.28p changes the style of interaction in the model. It depends on distance of oscillators 
to thresholds. We use this to introduce delay and continuous couplings in papers 42 1 and 43(, respectively. 

To illustrate Theorem 13.11 consider a group of oscillators, xi,i = 1,2,..., 100, with random uniform 
distributed start values in [0, 1]. It is supposed that they satisfy the equations x[ — (3+0.01/2^) — (2+0.01Ci)£i. 
The constants p,i,Ci, as well as £j in the thresholds 1 + 0.005^,i = 1,2,..., 100, are uniform random 
distributed numbers from [0,1]. In Figure [3] one can see the result of simulation with e = 0.08, where the 
state of the system is shown before the first, twenty first, forty second and sixty third firing of the system. So, 
it is obvious that eventually the model shows synchrony. Let us describe a more general system of oscillators 




20 40 60 BO 100 20 40 60 BO 100 



oscillators 



Figure 3: The state of the model before the first, twenty first, forty second and sixty third firing of the system. The flat sections 
of the graph are groups of synchronized oscillators. (Color online) 

such that Theorem 13. II is still true. A system of n oscillators is given, such that if z— th oscillator does not 
fire or jump up, then it satisfies the i— th equation of system If several oscillators Xi s , s — 1, 2, . . . , k, 

fire so that Xi a (t) — 1 + Q s (x) and x,i s (t+) = 0, then all other oscillators xi ,p = k + 1, k + 1, . . . , n, change 
their coordinates by the law 

x . (t+) = | 0, if x ip {t) + e + Eti <*,i. > 1 + dp{x), 
lp \ x ip (t) + e + J2 s =i e i P i s > otherwise. 

One can easily see that the last theorem is correct for the model just described, if e + 53 s =i e 'p i s > ^> ^ or au 
possible k, i p and i s . 

Remark 3.2. The analysis of non-identical oscillators with non-small parameters may shed light on the 
investigation of arrhythmias, chaotic flashing of fireflies, etc. Namely, the dynamics in the neighborhood 
and inside the periodic trajectory, Q in Figure [7J can be very complex. We do not exclude the possibility of 
chaos and fractals fl(\ l. Bifurcation of periodic solutions can be discussed, if the parameters are small. 

4. The Kamke condition and synchronization 

In this section we consider an integrate-and-fire model with a new type of continuous connection. 
We believe that models of identical oscillators with more general differential equations, 

x\ = g(xi,x 2 , ■ ■ -,x n ), 
x 'l = g[x2,x 3 , . . .,Xt), 



g(x n , X\ , . 
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(4.29) 



where < Xi < l,i = 1,2, ... ,n, are of both theoretical and applied interest. The positive valued func- 
tion g(y\, 2/2) • ■ • i Dn) hi (14.291) is continuously differentiable and indifferent with respect to permutations of 
coordinates y 2 to y n . 

When the oscillator Xj fires at the moment t such that Xj(t) = l,Xj(t+) = 0, then the value of an 
oscillator with i ^ j, changes so that 

*<«+>- {^—. (4 30) 

Consider the cone R™ C l n of all vectors with nonnegative coordinates. Introduce a partial order in the 
cone such that a < b if a, < 6,-, i = 1, 2, . . . , n, j46| . 

We say that function g is of type JC in R™ if a < b implies that g(a) < g(b). The sufficient condition for 

that is >0,i^l. 

Let u(t,t ,uo) and u(t, t , m), u , m £ R", be solutions of (|4.29p . If g is of type JC, then 
dynamics of (|4.29[) is monotone for t > 0. That is, u(t, 0, «o) < 0, «i), if u < «i. 

Consider first the model of two oscillators. Define map L for this system in the following way. Take the 
solution u(t) = u(t,0, (0,v + e)) = (1/1,1*2)- Denote by s(v) the moment when U2(s) — 1, and define the 
function L(v) = u%(s) on (0, 1 — e). Then define map L through (I2.6[) . Let us check if conditions (Al), (A3) 
are valid for this map. Indeed, the continuity of L is obvious. It is non-increasing since the monotonicity. 
Assume that there exist numbers vx, V2 € (0, 1 — e) such that v\ < V2 and s — s(vx) — s(v2)- Then, we have 
a contradiction as the open interval (vi,v 2 ) is mapped to the closed set {1}. That is, L satisfies {Al). 

Condition (A3) is easily verifiable. Now, one can determine sets Si, similarly to that in Section [5J and 
prove that the following theorem is valid. 

Theorem 4.1. Assume that g is of JC type, (A2) is valid, and t$ > is a firing moment such that Xi(to) = 
l,x±(to+) — 0. If x 2 (to+) € S m ,m is a natural number, then the couple X\,x-x synchronizes within the time 
interval [to + ^tT, to + mT\. 

We have, moreover, that, if g is of JC type and (A2) is not true then the system does not synchronize. 

Consider the multidimensional system of oscillators. Introduce the function G(y, z) = g(y, z,z, . . . , z), 
and define the integrate-and-fire model of two identical oscillators y and z with the following system of 
differential equations 

V 1 = G(y,z), 

z' = G(y,z). (4.31) 

Denote by u = (y, z), u(t) = u(t, 0, (0, v + e)), the solution of (|4.31[) . and by s(v) the moment when z(s) = 1. 
Next, define the function L(v) = y(s) on (0, 1 — e). Then map L can be defined by (|2.6p as well as correspond 
sets Si- By applying the monotonicity of the dynamics, one can prove the following assertion, in a way very 
similar to that of Theorem 13.11 

Theorem 4.2. Assume that < < n,i ^ 1, condition (A2) is valid, and to > is a firing moment 

such that Xj(to) — l,Xj(to+) = 0. If parameter 7/ is sufficiently small then the motion x(t) of the system 
synchronizes within the time interval [io,to + T], if Xi(to+) € So,* 7^ j, and within the time interval [to + 
max '^ 2 J — '- — -T, to + (maxi^j ki + 1)T], if there exist x s (to+) € [do, <2i) for some s ^ j and Xi(to+) G Sk, ,i 7^ j. 

Example 12.41 shows that the analysis of the map L for Theorems 14.11 and 14.21 is not simple even with linear 
differential equations. The results of this section are therefore are provided for numerical application, as 
well as for future investigations. 

5. Conclusion 

A version of the integrate-and-fire model of pulse-coupled and non-identical oscillators is investigated in 
this paper. We have made significant advances for the solution of second Peskin's conjecture, though we have 
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not showed that the measure of non-synchronized initial values is zero as it was shown in [2l| for identical 
oscillators. However, we have located non-synchronized points, and showed that the time of synchronization 
infinitely increases as the region of points to be synchronized enlarges. Moreover, it is not necessary, in 
applications, for all points of the domain to fire in unison, and it is sufficient for the neighborhood of a 
motion to be synchronized. 

A prototype map is introduced that helps to precise the results for a system of two identical oscillators and 
to solve the problem for the multi-non-identical-oscillators model. The approach of the paper is universal. 
For example, it can be applied if the thresholds are of the form 1 + cf>(t, x) or 1 + <f>(t), which represents an 
oscillating signal in physiology and a variation of the threshold of the electronic relaxation oscillator [32j • 
Moreover, the jump's value e as well as its perturbations may depend on x. The method can be easily extended 
for models, where differential equations have discontinuous right-hand sides, and only the existence of the 
map L with request properties are important to get appropriate results. Not only exhibitory, but inhibitory 
models of oscillators can be investigated, as well as delay of couplings 0, 0, 0- UK , and continuous couplings 
gen erated by firing (2BL 143^ . One can consider the models of our paper as discontinuous cooperative systems 
[46j |. Consequently, we expect that by applying methods of dynamical systems with variable moments of 



discontinuity [47| , the results on monotone systems [46|, |48|, |49| can be extended for these models 



There is a rich collection of results on oscillators, obtained through experiments and simulations. The 
approach of the present paper can give theoretical background for them and also form a basis for new ones. 
It can be applied not only to the problems of synchronization, but also to periodic, almost periodic motions, 
and the complex behavior of biological models. New small-world phenomena can be discovered. One can 
now request a certain property for the map L and then look for a system which meets the property. Thus, 
many new theoretical challenges can be brought under discussion. Conversely, if a system is given, then one 
can construct the corresponding map L, and by analyzing find new features which have not been mentioned 
in the present investigation. 

Acknowledgement 5.1. The author thanks Mehmet Turan and Zahid Samancioglu for the technical as- 
sistance. 
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